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ABSTRACT 


The oceanic turbulent boundary layer is a critical region to 
understand for oceanic and atmospheric prediction. This thesis 
answers two fundamental questions: (1) what is the response of the 
ocean mixed layer system to transient forcing at the air sea 
surface? (ii) what is the necessary time and space resolution in an 
ocean mixed layer model to resolve important transient responses? 

Beginning with replication of de Szoeke and Rhines’ (1976) 
work, additional physical processes were added to include more 
realistic viscous disSipation and anisotropy in the three- 
dimensional turbulent kinetic energy (TKE) budget. These refinements 
resulted in modification of de Szoeke and Rhines’ findings. 
Firstly, TKE unsteadiness is important for a minimum of 10° seconds. 
Secondly, viscous dissipation should not be approximated as simply 
proportional to shear production. Thirdly, entrainment shear 
production remains Significant for a minimum of one pendulum-day. 

The required temporal model resolution is dependent on the 
phenomena to be studied. This study focused on the diurnal, 
synoptic, and annual cycles, which the one-hour time step of the 
Naval Postgraduate School model adequately resolves. The study of 
Spatial resolution showed unexpectedly that model skill was 


comparable for 1m, 10 m and even 20 m vertical grid spacing. 
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I. INTRODUCTION 
A. STATEMENT OF PROBLEM 


The oceanic turbulent boundary layer or mixed layer is a 
critical region to understand when studying the atmosphere or 
ocean. The ocean is the primary source of moisture fluxed 
into the atmosphere, it is a significant source of heat fluxed 
into the atmosphere in many regions, and it is a moderator of 
atmospheric extremes. To accurately model the atmosphere or 
ocean, a realistic representation of the turbulent mixed layer 
is required. This thesis will address two fundamental 
questions: (i) what is the response of the ocean mixed layer 
system to transient forcing at the air sea surface?; (ii) what 
is the necessary time and space resolution in an ocean mixed 
layer model to resolve important transient responses? With 
the answers to these questions, air-ocean models will be a 
step closer to being realistically coupled. The result will 
be better predictions for both atmosphere and ocean. 

Be METHOD 

To accomplish this task, the work of de Szoeke and 
Rhines (1976) on wind driven mixing was replicated and then 
generalized to include surface heating and cooling. Their 
study in “Asymptotic regimes in mixed layer deepening” did not 
include physical processes associated with surface heat flux 


and anisotropy in the three-dimensional turbulent kinetic 


energy(TKE) budget. With the addition of these critical 
physical processes, a more generalized study of the short time 
scale dependence of the mixed layer is possible. The study 
can be extended past the single day scale of de Szoeke and 
Rhines to the synoptic and annual time scales by shifting the 
method of solution from a Runge-Kutta numerical scheme to a 
FORTRAN-based gridded model. With the guestion of the minimum 
time scale resolved, new assumptions about mixed layer physics 
can be made which allow testing of spatial resolution. The 
Naval Postgraduate School(NPS) mixed layer model developed by 
Garwood (1977) was modified to test the effect on model skill 
eu varying the vertical ied size using observed 
meteorological forcing and changing thermal structure at 


weather station Papa (50° N, 145° W). 


II. THEORY, MODEL EQUATIONS, AND HYPOTHETICAL SOLUTIONS 


A. AIR OCEAN COUPLING 


The surface layer of the ocean responds to variations in 
wind stress, heat flux, moisture flux, and several other less 
Pgiiiicanme —factors - These factors have been examined in 
depth since Ekman’s (1905)ground laying work. 

B. REPLICATION OF DE SZOEKE AND RHINES (1976) 

init@etals modeling “eliO@es, LOrm™enis study “Focused “on 
reproducing and then building on the efforts of de Szoeke and 
Rhines (1976) as discussed in their paper “Asymptotic regimes 
in mixed-layer deepening.” Niiler (1975) proposed the basic 
terms in the turbulent kinetic energy (TKE) equation studied 


by de Szoeke and Rhines: 


a 


(> NPht ul a 


S (Locos ft) +C,uch*) = =2myu:h (1) 


A B C D 
Term A represents buoyant damping of TKE by entrainment, term 
B the shear production of TKE by entrainment, term C the spin 
up or storage of TKE, and term D the near-surface wind shear 
PEOOUCtTION Of TRE minus Viscous dissipation. Solar heating 
and surface heat flux were not included. Term C, the energy 
required to spin up the turbulence intensity, was not 


eurginaliy ineludea by Niiler. 


To mathematically define the components of the de Szoeke 
and Rhines model, the friction velocity or u. must first be 
defined. The friction velocity (equation 2) is defined as 
the square root of the surface wind stress divided by the 


density of water. 





The Brunt-Vaisala or buoyancy frequency (N) is a measure of 
vertical stability. A value of 21/600 sec™’ was used for 
initial model testing and in the hypothetical cases. The 
mixed layer depth is h. The Coriolis parameter (f) is defined 
as 2Q, where Q, 1s the vertical component of planetary 
rotation. A nominal mid-latitude value of 1x107* sec” was 
used for f. In equation (1) C, and m are dimensionless tuning 
GOCELTCTENtS tor EXE @cpomneuoNanead nee wimnad-sShear produce ren, 
respectively. 

Buoyant damping of TKE by entrainment, term A, is the 
rate of decrease in turbulence due to the lifting upward of 
denser water originating below the mixed layer. The mixed 
layer is made more stable by the entrainment of colder denser 
water, with a commensurate increase in potential energy. 
Shear production due to entrainment, term B, is the conversion 


of mean kinetic energy to turbulent kinetic energy due to 


shear across the base of the mixed layer. This term also 
results in a net gain in turbulence. Turbulence spin up, term 
C, is the TKE build up required before entrainment can occur. 
The wind stirring minuS viscous dissipation, term OD, 
represents the net input of energy into the mixed layer due 
to wind stress, T,, in equation (2). In de Szoeke and Rhines’ 
model, viscous dissipation is also proportional to u.° because 
surface heating and cooling are neglected as sources or sinks 
of TKE, and dissipation of TKE produced by entrainment shear 
MPecocuction iS not Maeeunted ioe . 

As a first step in demonstrating and understanding the 
physics that drive mixed layer deepening, de Szoeke and 
Rhines’ results were replicated. Equation (1) was solved for 
the rate of mixed layer deepening, W,= oh/dt. The result is 
referred to as the W, equation, since the vertical change in 
mers due only to entrainment (no advection). The W. equation 
was formulated as a MATLAB function, and solved using an 
imercansic Runge-Kutta solution routine, ODE45. The MATLAB 
PMmeLllon iS restricted to continuous functions. This 
restricts the model to deepening scenarios. 

De Szoeke and Rhines simplified equation (1) by analyzing 
the time evolution of the mixed layer deepening process and 
determined which terms could be neglected during various time 


regimes. PWielvere Solutions to the equation wore msemon 


possible in those specific time regimes. Figure 2.la isa 
reproduction of de Szoeke and Rhines results, and 2.1b depicts 
ENeeresults Cf as nunge-huera scOoOlMuEron LO =-equation =). 

De Szoeke and Rhines showed that initially shear 
production due to entrainment (term B) can be ignored, that 
turbulence spin up (term C) dominates, and that the mixed 
layer deepens at a rate proportional to u.t. This condition 
holds for the first 100 seconds of deepening. Next, surface 
wind stirring (term D) dominates. Deepening continues at a 
Cates peopem=amonal to. t' =, Term C begins to build in 
importance between 20 and 120 minutes, with the deepening rate 
proportional to t’’*. Coriolis (term B) becomes a factor as the 
mean flow direction is turned away from the wind direction. 
After twelve hours, term B is no longer dominant. Deepening 
returns to a wind driven scenario and a rate proportional to 


ti/3 





Vo sec 


Figure 2.la de Szoeke and Rhines solution. Syeplaticl verb ig te © 


Hati1O D/A; dashed curve: ratio B/A; dotted curve: C/A 











time (seconds) 


pugube 2,15 MATLAB ODE45 solution, replicating de Szoeke and 


Rhines. 


ce IMPROVED MODEL PHYSICS TO INCLUDE EXPLICIT TKE 

In the second model, the basic W, equation is the same as 
in the first model, but three additional equations are 
included: the mixed layer temperature equation and the two 
horizontal mixed layer momentum equations. This allows the 
explicit calculation of the variables of temperature and 


velocity. The full set of equations are: 











= = . (e5") 
ale (Ee h 
CE i ACEO 
dt Oh h 
AVW 
ON 2 ys anne (5) 
dt Oh h 
es mass 
Aina Bien Ja 
ore 

= 1 2 : ae S 
C ee us — (1-cos ft) +Cyulh?) 


In this more complete model with explicit momentum and heat 
equations, the ability to add steady state upwelling was 
included in the W, equation. This model is identical to the 
Original de Szoeke and Rhines model if the upwelling term, 
W(z=-h), is set to zero, and if the net surface heat flux(Q,) 
is neglected. However, when Q, iS nonzero in the second 
model, Q, directly affects both the temperature equation (3) 


and the mixed layer depth equation (6), and it indirectly 


eee momo meliimncecalce Or the W. termS in equations (4) and 
eo eee eilomentumscquations, i, and t, are components of 
wind stress. AT is the change in temperature across the mixed 
Piet iectEaece Mi Piguges 2.2), Similarly, AU and,AV are the 


changes in the velocity components across the interface. 
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Megure 2.2 Hypothetical temperature profile. 

The second model’s solution for equation (6) was 
identical to de Szoeke and Rhines when Q and W(z=-h) were set 
to zero. The three additional equations allowed the explicit 
calculation of mixed layer temperature and velocity, but did 
not change the W, solution. The model was initialized with 


both horizontal velocity components equal to zero and an SST 


of 17° C with a linear temperature gradient with depth, having 
N = 21/600 seconds. A steady wind stress of 0.1 N/m* was 
applied in the 1, direction. In 30 hours, the SST decreased 
to 16.3 C (Figure 2.3) as the cooler sub-mixed layer water was 
entrained. The horizontal velocities, u and v, respond to the 
steady state wind stress in the form of forced inertial 
motion. The period is dependent on latitude (2n/f) but a 90 
degree phase shift between u and v is constant. 

The third model is a MATLAB version of the NPS mixed 
layer modell; it represents a Significant addition of TRE 
physics to the first two models. Equations for -atunmam 
heating and cooling and three for energy transformation were 
added, and the W, equation was reformulated. The resultant 


equations are: 











7 ee 

= Om iOum (ia x Oe eos ee 
Oe CL ICN 3000 ( oaaCay e 

Guth _» 9 FU? Wo 42m, (B-3u7) /E-= 25 FRE ie 
WE = 2, : m., ( a) es mM. ) 
Ovih yt +2 39°) (gee 2 ie fee 

a a Peis 3 E = (mE tide 5/073) 

oy ss 

Gw7h aghQ, —= 2 = 

Sa Ge) et eee — 3 Wy =— (nee +m CRE (10) 

a g EG, Bs ) VE = (my . ) 
No eee! eee (11) 
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A Q,(t) equation (7) was designed to have a diurnal structure 
with a net zero heat flux. The energy transformation terms 
Cones redistribute energy among components, representing 
pressure-strain interactions. The variable E is the summation 
of the squares of the turbulent velocity fluctuations in the 
<7. you anduez directions, F=u*+v+w . In equation (8), mete 
is the near surface wind shear production term. The summation 
of the three m,E°’? + mfhE terms found in equations (8), (9), 
and (10) is the viscous dissipation. The m,u.* term and the 
summed m,B°/? + mfhE terms are approximated in the de Szoeke 
and Rhines equation by the mu.-h term. The m,(E-3u 7) E/” terms 
are pressure redistribution of TKE among the components, in 
which the turbulence tends toward isotropy by pressure-strain 
interactions. The aghATW, term in equation (10) is the 
buoyant damping of TKE by entrainment and corresponds to de 
Szoeke and Rhines N*h* term. The aghQ,/pC, term in equation 
(10) is the buoyant damping of TKE by surface heating (Q,), 
and has no corresponding term in equation (1), although it was 
included in equation (6) of the second model. Equation (11) 
is the entrainment rate as a function of the ambient TKE. 
Equation (12) is the rate of change of the mixed layer due to 


entrainment, W., and vertical mean advection by W(z=-h). 


TZ 


Initial runs of this model were conducted with a constant 
wind stress but no heating or cooling. These solutions were 
Similiar to the results of models one and two, but the mixed 
iPjer Cepens Were ySigniticantly shallower (Figure 2.4) > has 
is as expected and is attributed to energy being expended for 
turbulence spin up and dissipation in the layer as well as 
entrainment. Two differences between the first and third 
models are significant: (i) calculation of dissipation and 
Ge) role Of iImertial oscillations. 

Term D, in the de Szoeke and Rhines model, near~surface 
wind shear production of TKE minus viscous dissipation, is an 
ever Simplaetication. ihe viscous dissipation sheulld net be 
represented as directly proportional to u.*h. The dissipation 
resulting from the shear production due to entrainment is not 
represented at all. Figure 2.5 is analogous Figure 2.1, and 
me shows the spin up Of model terms. The ratio of O03/A3 is 
mearc-surface wind shear production of TKE minus viseous 
dissipation over buoyant damping of TKE by entrainment. The 
ratio of B3/A3 is shear production of TKE by entrainment over 
buoyant damping of TKE by entrainment. The ratio C3/A3 is the 
turbulence spin-up term over buoyant damping of TKE by 
entrainment. The most obvious difference is in the ratio 
D3/A3. Total dissipation is much larger in the third model 


than with the approximation used in the first model because 
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Figure 2.4 Comparison of mixed layer depths between model two 


and model three. 
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Pirgure. 2.2 Model three solution. The ratios represent the 


same physical processes as Figures 2.1a and agls O@ 


dies, 


ENO first Model Mest imared VISCOUS GUISsSsipacion as) preperer oma 
to the wind stress only. The third model calculates the 
dissipation for all three TKE components based on the 
explicitly calculated TKE in the mixed layer. Between 10° 
seconds and 4x10* seconds, there is more viscous dissipation 
than wind shear production, because of the addition of the 
etissipatten due to shear production: 

The significance of term B decreases as the mixed layer 
deepens. In the first model, the effect of inertial forcing 
is negligible after 6x10* seconds. In the third model, term 
B3 remains significant four times longer than for model one. 
If the dissipation due to shear entrainment is neglected, the 
amplitude of term B3 decreases, and term A3 and D3 dominate 
(Figure 2.6) after 1x10° seconds. 

D. MODEL RESOLUTION IN THE TIME DOMAIN 

Models one through three all use a Runge-Kutta solution 
with a variable time step to solve the mixed layer equations. 
The numerical scheme maximizes the time step to solve the set 
of equations without exceeding specified error tolerances. [In 
mixed layer physics, oh/dt is not a continuous function when 
mixed layer shallowing occurs. For this reason, the Runge- 
Kutta solution can only be used during periods of mixed layer 
deepening. To study realistic cases, a model that can shallow 


as well as deepen was required. 


iG 
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Bigure 2.6 Model three solution wath viscous dissipation 


Coriolis parameter m,=0. 


ie 


The NPS mixed layer model was chosen. The NPS model is 
a FORTRAN-based gridded numerical model with a fixed time step 
of one hour. Using a fixed time step to simultaneously solve 
differential equations introduces errors. To determine the 
adequacy of the one-hour time step, the NPS model was compared 
to the third model using mixed layer depth (Figure 2.7) and 
temperature (Figure 2.8) as the basis of comparison. The NPS 
model does not have a term for shear production of TKE by 
entrainment (term B3) but vertical mixing by dynamze 
instability occurs when the Richardson number (Ri) is less 
Eneameewe Critical value, Ri < 0.25, following Che procedure op 
Adamec et al. (1981). These physical processes can then be 
selectively included or neglected for comparison purposes. 

The NPS model deepens too gquickly during the first time 
step. This 1s because the NPS model assumes fully spun-up 
steady state turbulence. This is accomplished by solving the 
steady state equations (13), (14), and (15), or by setting 


time lett hand side of equations (8), (9), and (10) to zero. 


0=V°W,+2m, (E- 3 v2) yE-< (x 2m. £ hE) les) 
aghQ —_ ie 
0 =-aghA TW,-——* +2 m, (E-3 w*) YE- = (m, E ? +m, £hE) Ce) 
O 
p 
2 
0=2m,u;+U?W+2m, ie 3u? ) YE- = “+m. £hE) (15) 
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Pogure 2.7 Mixed layer depth comparison between the NPS mixed 
layer model with Richardson Number deepening turned off and 


mecel mumaer senree with Term BS turned OF. 
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- — Model 3 





’ 
76.1 16.2 16.3 16.4 Wee) 16.6 1657 16.8 


Temperature (degrees C) 


Figure 2.8 Mixed layer temperature comparison between the NPS 


mixed layer model and model number three, at time = 8 hours. 


Z10 


This assumption allows instantaneous deepening to occur 
without having to spin up the turbulence. Depending on the 
faecmOnerlNCEUg@elGn im Ehe forcing, this may not be a valid 
Sescupe tone 2toGure Z2Zaemgepicks a portion of the C3/A3 ratio 
‘SieGiel  Sak(e @bhat= e742) Examination of the results depicted in 
Figure 2.9 reveals that the unsteady term does not approach 
zero at 10* seconds as proposed by de Szoeke and Rhines but at 
meeescconds. At 0.36x%10* seconds, this assumption results in 
a maximum 14 percent error in TKE generation for any change in 
Bemcing. If the Eorcing 1s constant or slowly changang this 
error is negligible. Figure 2.9 indicates that the unsteady 
term persists and oscillates in excess of 3 days. 

The temperature profiles differ for two reasons. The 
first is related to mixed layer depth. The NPS mixed layer 
model depth is shallower and therefore warmer. Conversely, 
model three has a deeper mixed layer depth and is cooler. The 
second reason is related to model thermal structure. Model 
three has a continuous function for temperature versus depth. 
The NPS model is a gridded model and assumes a step like 
mMCELON LOr temperature. The @aifference in thermal structure 
between the models is evident in Figure 2.8 in the region 
below the mixed layer. 

To further test the adequacy of the one hour time step, 


the NPS model was run with a 15 minute time step. The results 
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Figure 2.9 Ratio of Model three turbulence spin-up over 


buoyant damping of TKE by entrainment (C3/A3). 


ae 


are depicted in Figure 2.10. After the second time step of 
the one hour model, the mixed layer depths of the two NPS 
models are parallel in time with an offset of approximately 
0.5 m. The one-hour time step of the NPS mixed layer model is 
sufficient for the proposed study. This conclusion was 
reached because the error induced by the relatively large time 
step is small and explainable. After initial spin up, the 15 
minute time step does not result in significant resolution 
improvement. Additionally, the subsequent forcing data to be 
used has a sampling period of three hours. Interpolating to 
time steps less than an hour was therefore judged to be 


inappropriate. 
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Figure 2.10 Mixed layer depth comparison using the one-hour 
time step and the 15 minute time step NPS mixed layer model, 


and model number three. 
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III. REALISTIC SOLUTIONS 


Doc IMPORTANCE OF VERTICAL RESOLUTION 

In any geophysical process model, the resolution required 
is dependent on the scale of the physical phenomena to be 
represented. The annual fluctuation in the mixed layer depth 
Sere ceatlom Papa “rs on the order Vor 150° m- The daily 
fluctuations vary between order 10 m and order 100 i, 
depending on time of year and the synoptic forcing. Unlike 
other mixed layer models, the NPS mixed layer model used here 
employs a floating grid point assigned to the mixed layer 
depth to increase the model’s resolution. This feature 
becomes more significant as grid size increases. 


The choice of vertical grid spacing (Az) is governed by 


several factors. The two factors considered in this thesis 
are model computational speed and resolution. These two 
factors are not independent of each other. Smaller grid 


spacing requires increased computer time, but resolves more 
features; larger grid spacing is faster, but has reduced 
resolution. Model requirements will determine which one of 
these two factors is more important. In the past, embedding 
a mixed layer within an ocean general circulation model or 
coupling a well-resolved upper ocean model to a global 


atmospheric model has been prohibitive because of the number 
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of grid levels required to represent the ocean. The following 
sections address the question of required vertical grid 
Spacing and the accuracy of larger grid sizes. 
B. MODEL FORCING 

The atmospheric forcing utilized was calculated from 
observations by a Canadian weather ship in the vicinity of 
50° North latitude 145° West longitude, nominally referred to 
as Station Papa. Three hourly weather observations were 
taken, consisting of wind speed and direction, air and sea 
Surface temperature, dew point, and fractional cloud cover. 
The three-hourly observations were interpolated to hourly 
values. Aerodynamic bulk formulas derived by Large and Pond 
(1982) and modified by Oberhuber (1993) were used to compute 
the momentum, sensible, and latent heat fluxes, correcting for 
atmospheric stability. Mechanical bathythermographs (MBT) were 
also taken at irregular intervals at station Papa, varying 
from less than an hour to several days. The MBT data was 
manually digitized at 5m intervals. The MBT data are used 
in this study to initialize the mixed layer model and then for 
subsequent verification. Model runs were conducted for 
various periods between 1961 and 1969. In addition to being 
one of the longest running continuous oceanic time series, the 


data set is ideally suited to testing a one dimensional model 
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because the horizontal currents are weak and advection due to 
upwelling is generally believed insignificant (Tabata, 1965). 
C. EFFECT OF VERTICAL GRID SIZE 

The 1m grid NPS mixed layer model has a minimum mixed 
imver depth of slammior sone Az). Changing the grid spacing 
raised the issue of what the minimum mixed layer depth or 
ceiling should be. When Az was increased to 10m and 20m, 
a 1/4 Az minimum mixed layer depth was found to be optimal. 
The 1/4 Az ceiling resulted in a 30% decrease in mean error in 
temperature over the one-Az ceiling. In the presence of 
strong heating and no wind stress, ceilings smaller than 1/4 
Az resulted in the mixed layer temperature becoming 
unrealistically hot. The same unrealistic temperatures were 
also obtained in the 1m grid case by raising the ceiling 
above one Az. 
D. VERIFICATION WITH OBSERVATIONS 

a. Model Predicted SST vs Bucket Temperature 

Sea surface temperature is the only state variable for 
which model output and station Papa observations can be 
@emeanucusly G€ompared: To accomplish this task, thew mode! 
mixed layer temperatures were compared to the interpolated 
bucket temperatures. The mean difference and standard 


deviation were used for comparison. 


Za 


ap ( ae > SST cket ) 
18) 


The mean difference = 


The mean difference calculation was used to evaluate the 
effect of varying the minimum allowable mixed layer depth as 
described in section III C. The results presented in Appendix 
A provide the justification for selecting 1/4 Az as the 
iin cel Ling Lom Az so Of 10 mand 20m. lhe elgme -anmuUalns 
period model runs commenced in mid-March, near the time of 
maximum mixed layer depth, and ran for 365 days. The results 
of the mean differences were unexpected and raised questions 
about optimal spatial resolution. In every case but one, the 
error for model runs with 1m grid spacing was greater than 
for both the 10 m and 20 m Az cases. This implies that the 
coarser grid is more accurate. The standard deviation showed 
the same trend. The model runs depicted in Figure 3.1 are 
characteristic of the other years. Spectra were calculated 
for lem, lLO0=m encwreeUenegrid resolution (Proqures 32,333, 4 
3.4, respectively). The annual cycle was not resolved due to 
the spectral computational length of the time series; however, 
the diurnal cycle is evident. The diurnal cycle was best 
resolved by the 10 m grid, then by the 20 m grid, and lastly 
Dyeene 1 mm Grid. 

In an attempt to understand why better model performance 


is achieved with coarse resolution, model tuning sensitivity 
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tests were conducted. Tuning coefficients P3 and AM3 in 
equations (13), (14), and (15) were varied, and the model 
Zo LiierOorseichwyear., Tuning coefficient P3 equals m./m,, and 
feGuUudueomencmerteockh Or the Coriolis in the dissipation term; 
itagen PS 9reSumes in Stronger dissipation. Tuning coctficienr 
AM3 equals m; and regulates the efficiency of shear 
production. Larger AM3 results in increased shear production. 

Year-long model runs for eight years were made for each 
P3 and AM3 pair. The values ranged for 3 to 8 for P3 and 1 to 
9 for AM3. The standard deviation for each run was calculated 
and the results averaged. Figure 3.5 depicts the average 
eaceor ror the | m grid spacing cases; Figures 3.6 and 3.7 
depict the average error for the 10 m and 20 m grid spacing 
respectively. A well-defined region of minimum error is 
evident for the Az equal to 10 m and 20 m cases. This same 
region.is not as clearly defined in the Az equal to 1m caSe. 

Optimal tuning parameters were chosen which represent the 
eo, AMS) region of minimum error. The model was rerun for 
five different time frames and for seven different tuning 
coefficient pairs. The first was a repeat of the annual 
cycle. The next four were 10 day runs, one per season. This 
was done not only to elucidate differences in model accuracy 
among the seasons, but also to determine if the 1 m Az model 


improved over the larger grid models for diurnal and synoptic 
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scales. The standard deviations were calculated and 
tabulated. Appendix B is the complete listing of results by 
tuntmg factor pairs. 

The results of the in-depth tuning study are both 
Seorlevying and conclusive. the vertical grid spacing accumadey 
is dependent on time of year and length of model run (see 
table 1). During the spring and summer when the forcing is 
dominated by the diurnal cycle, the larger Az’s are more 
accurate. During the synoptically forced deepening period 
of the fall, the 1m grid spacing more accurately represents 
the bucket temperatures. In the 48 cases where the grid size 
is a factor when using 10 m vice 1m Az, an increase in 
standard deviation error of 0.06° C resulted, and an increase 
of 0.05° C resulted in the 35 cases when 20 m Az were used. 
The winter accuracy is generally independent of grid spacing. 

Z: Bucket vs Mechanical Bathythermograph SST 

The NPS mixed layer model uses BT’s to initialize the 
model . After Inverell Zzat Lom, surface meteorological 
observations are used to drive the model. The surface 
observations were taken at 3 hourly intervals, but the MBTs 
were deployed at irregular intervals. For initial conditions, 
the model uses the BT with the closest Julian date to the 
start date of the run. In general the relative vertical 


structure of the surface layer will be correct, but a random 


Sy 


bias in the absolute temperature will be induced 
Additionally, a phase shift is also likely, as the model 
starts at 00 GMT and the BT may or may not be at QO GMT. 
Figure 3.8 demonstrates the bias between bucket and BT 


temperatures and the irregular sampling interval of the BTs. 
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IV. SUMMARY AND RECOMMENDATIONS 


This study has evaluated the response of the ocean mixed 
layer to transient forcing and the time and space resolution 
aie rOucdceccuracely model upper ocean Variability. ahis 
study began with replication of de Szoeke and Rhines (1976) 
work, using a Runge-Kutta MATLAB solution method. Additional 
physics were added to their model including: surface heat 
flux, improved viscous dissipation, and anisotropy in the 
three-dimensional turbulent kinetic energy (TKE) budget. The 
surface heat flux was added for completeness but was not 
extensively tested in evaluating the most rapid transient 
responses with a less than diurnal period. The refinement of 
Viscous dissipation and TKE budget resulted in modification of 
de Szoeke and Rhines’ findings. First, the unsteady TKE term 
is important for a minimum of 10° seconds rather than only 10% 
seconds. Next, viscous dissipation should not be approximated 
as being proportional simply to u.*h. This parameterization 
fails to take into account buoyancy flux due to surface heat 
flux and entrainment. The time scale of the turbulence also 
@lays a role in dissipation, and 1S not previously included. 
Finally, shear production due to entrainment remains 
Significant for a minimum of one pendulum-day instead of the 


one-half pendulum-day of the de Szoeke and Rhines’ study. 
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Expanding the investigation to model predictions longer 
than a day, the NPS FORTRAN mixed layer model was employed. 
The required temporal resolution of mixed layer models is 
dependent on the phenomena to be studied. This study focused 
on the diurnal, synoptic, and annual cycles, for which the 
one-hour time step of the NPS FORTRAN model is more than 
adequate to resolve. The one-hour time step is well matched 
with the larger vertical grid spacings. Further research is 
required to determine whether a one-hour time step is too 
coarse when used with 1m grid spacing. If a time step 
shorter than one hour is to be used, the unsteady term may 
need to be added to the NPS FORTRAN model. 

Finadaly, the question of spatial resolution was 
addressed. The initial hypothesis for this research was that 
smaller grid spacing would yield better results than larger 
grid spacing. The results reported in Appendix A and Table 1 
were therefore unexpected. With the exception of the winter 
time frame, the 10 m and 20 m grid spacing generally equaled 
OG OMt—perrormed the Tem grid ssacing-, Use of the 20 m grid 
Spacing during Eneewinter resulted imgan rmsSHenuror of 0.05° CG, 
which is smaller than the Q0.1° C precision of the bucket 
temperatures. Overall, the NPS FORTRAN model skill was 
comparable for 1m, 10 m, and even 20 m vertical grid spacing. 


Therefore, the conclusion of this thesis is that use of 10 m 
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dicmezO Mevercercalegqrid Spacing in global ocean Models and 
coupled air-ocean models to represent the mixed layer is 


accurate and justified. 
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APPENDIX A. ANNUAL MEAN OF PREDICTED TEMPERATURE ERROR (DEGREES C) FOR 
VARIOUS MINIMUM ALLOWABLE MIXED LAYER DEPTH 
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APPENDIX B. RMS TEMPERATURE ERROR (DEGREES C) VERSUS MODEL AZ AND 
TUNNING PARAMETERS P3 AND AM3 
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